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Abstract 

We investigate the possible influence of the MSW effect on the expectations 
for the solar neutrino experiments in the maximal mixing scenario suggested by 
the atmospheric neutrino data. A direct numerical calculation of matter induced 
effects in the Sun shows that the naive vacuum predictions are left completely 
undisturbed in the particular case of maximal mixing, so that the MSW effect 
turns out to be unobservable. We give a qualitative explanation of this result. 
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There is no doubt that the famous Mikheyev-Smirnov-Wolfenstein (MSW) mecha- 
nism Hlj] continues to provide an elegant and viable explanation for the existing 
solar neutrino data [§] [0- The preferred MSW fit requires one neutrino mass- 
squared difference ~ 10~^ eV^ and one small mixing angle ~ 10~^ radians. On the 
other hand large mixing, and in particular maximal mixing [Q is not completely ruled 
out by the existing solar data and is in fact actively suggested by independent data 
relating to atmospheric neutrinos |10| [|ri|] ||12|. The atmospheric neutrino data 
require a larger neutrino mass-squared difference ~ 10~^ eV^. In this paper we focus 



attention on the maximal mixing scenario |J13[ ||14| suggested by the atmospheric data 
and proceed to investigate the possible influence of the MSW effect on the expecta- 
tions for the solar neutrino experiments in that case. Our main results are based on a 
direct numerical calculation of matter induced effects in the Sun. We find that, in the 
maximal mixing scenario, matter induced effects turn out to be essentially unobserv- 
able, with the naive vacuum predictions left completely undisturbed, in the specific 
case of maximal mixing. 

The MSW effect has its origin in the interaction of the solar neutrinos with the 
matter in the Sun. In particular in the presence of matter the neutrino mass matrix is 
modified by the forward scattering of electron-neutrinos from electrons via the weak 
charged current. In a weak basis which diagonalises the mass matrix of the charged 
leptons, the 3x3 vacuum propagation matrix m? /2E is replaced by a matter propaga- 
tion matrix m^/2ii^-|-diag(v^GAi'e, 0, 0), where mm) = m^ is the hermitian-square of 
the vacuum neutrino mass matrix, E is the neutrino energy, G is the Fermi constant, 
and Nf. is the (position dependent) number density of electrons in the Sun. 

We calculate the MSW effect numerically for arbitrary 3x3 vacuum mixing. To 
specify the vacuum mixing we take over the standard parameterisation [|15| of the quark 
mixing matrix, so that, for example, threefold maximal mixing |13[ is reproduced by 
setting 012 = 623 = 7r/4, ^13 = sin"^(l/v/3) and 5 = n/2. For given values of the input 
mixing parameters, we first construct a vacuum neutrino mass matrix (in the above 
basis) as a function of the two independent neutrino mass-squared differences Am^ and 
Am'^ (Am^ > Am'^ ||13[). To account for matter effects we divide the propagation path 
longitudinally into thin slices of (variable) thickness A. For a given slice we calculate 
the matrix of transition amplitudes just as in the vacuum case {A = exp{—im'^A/2E)), 
but using the matter propagation matrix calculated assuming a constant density over 
the slice. The overall matrix of transition amplitudes is given by the ordered product 
of those for the individual slices, and the final electron-neutrino survival probability 
P(e e) is averaged over neutrino energies for comparison with experiment. 

Our calculation is 'exact' (at least in the limit A — > 0) in the sense that it 
does not rely on any particular physical approximation relating to the MSW effect 



itself (eg. 'adiabatic approximation' etc. ). In practice, for the results presented be- 
low, we consider only radially directed propagation paths starting from the center of 
the Sun (the detailed production profile is unimportant here) and we average over 
a uniform distribution of neutrino energies with a width of ±25% The num- 



ber density of electrons in the Sun as a function of radius R is parameterised [16 



by = Ar^exp(5.50 - ^{IOMR/Rq)^ + 0.892), ^^ere Na = 6.02 x lO^^ cm-^ and 
Rq = 0.7 X 10^ km is the solar radius. For given Am^ and Am'^, typically ~ 1000 
propagation slices and ~ 2500 energy samplings proved sufficient to produce a robust 
result. High precision (128-bit) arithmetic was found to be neccessary for carrying out 
the matrix manipulations. 

We first present results for the case of 2 x 2 mixing, ie. for a simplified scenario in 
which one of the three generations is completely decoupled in the mixing matrix and 
may in effect be forgotten. In the 2x2 case the mixing is completely specified by a 
single mass-squared difference Am^, and by a single mixing angle 6. Figure 1 shows 
our results for 2x2 mixing. Figure la is for vacuum mixing only, ie. with matter 
induced effects neglected, and Figure lb shows our results with matter induced effects 
taken into account. We plot the expected electron survival probability P(e e) for 
solar neutrinos, as measured on Earth, as a function of Arn^/E, for various values 
of sin^ as indicated. For Am'^/E ^ 10"^^ eV^/MeV, the oscillation length is longer 
than the distance from the Sun to the Earth and P{e — e) = 1. In the case of 
vacuum mixing the biggest suppression (a factor of 1/2) occurs in the case of twofold 
maximal mixing (sin 9 = l/\/2). With matter effects included, for sin 9 < 1/V2 (eg. for 
sin 6* = 0.5, see Figure lb) we reproduce the familiar MSW 'bathtub' suppression, 
extending over the range Am? / E = 10^^ — 10^^ eV^/MeV. For sin 9 > l/y/2 (eg. for 
sin^^ = 0.9, see Figure lb) we have an inverted bathtub, ie. an MSW enhancement. 
The MSW enhancement occurs when the lighter charged lepton, the electron, couples 
preferentially to the heavier neutrino (the rows and columns of the mixing matrix are 
ordered in increasing mass). In the case of twofold maximal mixing there is neither 
a suppression nor an enhancement. The solid curve in Figure lb is identical to the 
corresponding curve in Figure la. We conclude that in the 2x2 context the MSW 
effect is unobservable in the particular case of maximal mixing. 

Our results for 3x3 mixing are shown in Figure 2. For the 3x3 case we compute 
the survival probability P{e ^ e) as a function of the smaller independent mass- 
squared difference Am'^, with the larger mass-squared difference Am^ fixed by the 
atmospheric neutrino data (Am^ = 10"^ eV^). The MSW effect is now governed by 
\U2e\, the magnitude of the element of the mixing matrix linking the electron with 



the second lightest neutrino mass eigenstate. In the standard parameterisation [15 



U2e = cos sin 6^12. For the results presented in Figure 2 we vary \U2e\ by varying 9i2, 



keeping 6*23, 6*13 and 6 fixed. For Am''^/E < lO'^^ eVVMeV we have P(e e) = 5/9. 
Figure 2a is for vacuum mixing only and Figure 2b shows our results with matter 
effects included. In the case of vacuum mixing the biggest suppression (a factor of 
1/3 for 3x3 mixing) occurs in the case of threefold maximal mixing. With matter 
effects included, the MSW effect leads to a suppression or an enhancement, in general, 
depending on the value of ^12, as shown by the broken curves in Figure 2b. The MSW 
effect has no influence at all, however, in the particular case of threefold maximal 
mixing (cf. the solid curves in Figure 2b and Figure 2a) mirroring exactly our results 
for the 2x2 case above. 

For completeness it should be said that if the larger mass-squared difference Am^ 
were not constrained by the atmospheric data and if it fell in the appropriate range 
viz. Am^ = 10~^ — 10~^ eV^, then the MSW effect would become observable in the 
case of threefold maximal mixing. The effect, however (assuming Am^ ^ Am'^), is 
simply to suppress P(e — e) by a factor of 1/3 (instead of 5/9) over the range of the 
bathtub, as shown in Figure 3, so that again the observable suppression factors are in 
general identical to the case of vacimm mixing. 

With a view to obtaining a better understanding of the evident 'special case' sta- 
tus of maximal mixing with respect to the MSW effect, we return, for simplicity, to 
reconsider the case of 2 x 2 mixing, as a function of Am'^/E, for an arbitrary mixing 
angle 6, as before. The effective Hamiltonian for the MSW system is just the matter 
propagation matrix above, which in the 2x2 case — (say), may be written: 



In vacuum, in the small 9 limit, is the heavy mass eigenstate. In the high density 
hmit Ve is the heavy mass eigenstate. The famihar near-total MSW suppression of 
the Ve flux for small mixing angles occurs when the matter density proflle in the Sun 
provides a smooth matching from the state at the point of production, to a near-i/^ 
state outside the Sun. 

We exploit an analogy between the v^, — system in the presence of a variable 
matter density, and the behaviour of a spin-1/2 dipole at rest in a time-dependent 
uniform magnetic field. Suppose that the dipole has a negative magnetic moment — /x. 
Suppose further that the magnetic field B seen by the dipole may be decomposed as 
the vector sum of a constant (ie. time independent) 'intrinsic' field B*^ and a variable 
(ie. time dependent) 'external' field Bg. If the external field Bg is directed along 
the quantisation axis (the 2;-axis) while the intrinsic field B° makes an angle 29 with 
respect to the negative z-axis and is contained in the zx-plane, then the Hamiltonian 




(1) 



for the dipole may be written: 

/ cos 2^ + /iSe ^S0sin2^ \ , . 

V piB'^ sm2e iiB^ cos2e - iiB J ' 

Comparing Eqs. 1 and 2 we see that, with the correspondence /iB^ ^ Am^/(2£') 
and iiBe <-> GNe/V^, the dipole and the MSW system have the same Hamiltonian. 
The utihty of the analogy hes in the fact that the behaviour of the dipole is readily 
understood, since the spin vector S for the dipole satisfies a well known classical equa- 
tion of motion (S= — 2yuSAB). The dipole simply precesses around the instantaneous 
magnetic field B, with instantaneous angular frequency 2fiB. 

The small angle MSW effect may now be viewed as the adiabatic reversal of the 
spin of the dipole, in response to the slow reversal of the field B, as the external 
field Be decreases to zero. This is illustrated in Figure 4. The resonance condition 
is satisfied when B is directed horizontally along the x-axis, and the mixing becomes 
momentarily maximal. If the vacuum mixing is anyway maximal the intrinsic field B° 
is directed entirely horizontally along the a:-axis and no such reversal can occur. The 
inverted bathtub seen for large mixing angles corresponds to the case that B° points 
upwards. In the maximal mixing case, the dipole simply follows B smoothly from the 
vertical to the horizontal and remains horizontal, yielding 50% spin-up (z/g) and 50% 
spin-down (i/^u), just as for maximal mixing in vacuum. The difference is that in the 
matter case the residual oscillations are small, while in the vacuum case the dipole 
precesses around the x-axis, corresponding to maximal (100%) oscillations. It is only 
because these oscillations are unresolved that the matter and the vacuum predictions 
turn out to be indistinguishable. Presumably the behaviour in the 3x3 case has some 
closely related explanation, but we have not attempted to consider the 3x3 case in 
the equivalent level of detail. 

To summarise, a direct numerical calculation shows that the MSW effect is unob- 
servable in the particular case of maximal mixing. The naive vacuum predictions are 
left completely undisturbed in that case. This result is valid for 2x2 mixing for any 
value of Am^, and for 3x3 mixing for any value of Am'^, with Am^ ~ lO^^cV^, fixed 
by the atmospheric neutrino data. Exploiting the analogy between the MSW effect 
in 2 X 2 mixing and the behaviour of a spin- 1/2 dipole in a time-dependent magnetic 
field, we have given a qualitative explanation of this result. It is true nonetheless 
that the small angle MSW solution, with appropriate choices for the parameters, gives 
an excellent fit to the existing solar data, and is currently (perhaps not unnaturally 
therefore) widely accepted as the solution to the solar neutrino problem. While our 
results do not in any way undermine the validity of the MSW solution, they do serve 
to draw attention to an interesting and significant exception, where the MSW mecha- 
nism cannot be invoked. The MSW and maximal mixing scenarios as solutions to the 



solar neutrino problem are, in a very definite and real sense, to be seen as mutually 
exclusive alternatives. 
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Figure Captions 



Figure 1. The expected electron-neutrino survival probability P{e — > e), for solar 
neutrinos, as measured on Earth, for the case of 2 x 2 mixing. The survival probability 
is plotted as a function of Am^/ E for various values of sin 9, as indicated, for a) vacuum 
mixing and b) accounting for matter induced effects in the Sun. In the particular case 
of maximal mixing (solid curves) the vacuum and matter curves are indistinguishable. 

Figure 2. The expected electron-neutrino survival probability P(e — > e), for solar 
neutrinos, as measured on Earth, for the case of 3 x 3 mixing. The survival probability 
is plotted as a function of Am''^/E (with Am^ = 10~^ eV^) for various values of ^12, 
as indicated, for a) vacuum mixing and b) accounting for matter induced effects in 
the Sun. In the particular case of threefold maximal mixing (solid curves) the vacuum 
and matter curves are indistinguishable. 

Figure 3. If Am^ were not constrained by the atmospheric data, the MSW effect would 
become observable in threefold maximal mixing for Am^ = 10~^ — 10^^ eV^, as shown 
(assuming Am'^ <^ Am^). The suppression factor is 1/3 (instead of 5/9) over the 
range of the 'bathtub', however, so that the observable suppression factors are anyway 
identical to those for the vacuum case. 

Figure 4. The near-total MSW suppression of the Ue flux for small mixing angles is 
analogous to the adiabatic reversal of a spin- 1/2 dipole in the time dependent magnetic 
field illustrated. As the 'external' field Be decreases to zero for fixed 'intrinsic' field 
B° as shown, the resultant field B seen by the dipole reverses. In the case of maximal 
mixing the intrinsic field B° is directed entirely along the x-axis and no such reversal 
occurs. 
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